Interface dynamics in Hele-Shaw flows with centrifugal forces. 
Preventing cusp singularities with rotation 
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A class of exact solutions of Hele-Shaw flows without surface tension in a rotating cell is reported. 
We show that the interplay between injection and rotation modifies drastically the scenario of 
formation of finite-time cusp singularities. For a subclass of solutions, we show that, for any given 
initial condition, there exists a critical rotation rate above which cusp formation is prevented. We 
also find an exact sufficient condition to avoid cusps simultaneously for all initial conditions. This 
condition admits a simple interpretation related to the linear stability problem. 
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The dynamics of the interface between viscous fluids 
confined in a Hele-Shaw cell Q-pj has received attention 
for several decades from physicists, mathematicians and 
engineers. In particular it has played a central role in the 
context of interfacial pattern formation HJH]. As a free 
boundary problem it has the particular interest that ex- 
plicit time-dependent solutions can often be found in the 
case with no surface tension [§|-^. As a consequence, the 
issue of the role of surface tension as a singular pertur- 
bation in the interface dynamics has received increasing 
attention [ll^-H for its potential relevance to a broad 
class of problems. However, to what extent the physics 
of the real problem (with finite surface tension) is cap- 
tured, even at a qualitative level, by the known solutions 
is still poorly understood. 

As an initial-value problem, the zero surface tension 
case is known to be ill-posed [O. An important aspect 
related to this fact is that some smooth initial conditions 
develop finite-tinre singularities in the form of cusps of 
the interface &n|. After this blow-up of the solution, 
the time evolution is no longer defined. Generation of 
flnite-time singularities is in itself interesting in connec- 
tion with other singular perturbation problems in fluid 
dynamics, such as in the case of the Euler equations. In 
the present problem, surface tension acts as the natural 
regulator curing this singular behavior, but unfortunately 
the problem with surface tension is much more difflcult 
and usually defies the analytical treatment. Motivated by 
this fact, and inspired by recent experiments on rotating 
Hele-Shaw cells ||l^,|l^, we address here the 'perturba- 
tion' of the original free boundary problem by the pres- 
ence of a centrifugal field. This new ingredient enriches 
the problem in a nontrivial way but, as we will see, it 
still admits explicit solutions, and may thus lead to new 
analytical insights, both in understanding the generation 
of finite-time singularities, and in elucidating the dynam- 
ical mechanisms of Laplacian growth with and without 
surface tension. In fact, although the new dimensionless 
parameter introduced by rotation does not fully regular- 
ize the problem, we show that it may prevent the oc- 



currence of finite-time singularities, allowing for smooth, 
nonsingular solutions for the entire time evolution. This 
fact alone enlarges the class of (nontrivial) exact solutions 
without surface tension which are potentially relevant to 
the physically realizable situations. 

We study an interface between a fluid with viscosity n 
and density p and one with zero viscosity and zero den- 
sity in a Hele-Shaw cell with gap h. The cell can be put in 
rotation with angular velocity il and fluid can be injected 
or sucked out of the cell through a hole at the center of 
rotation, with areal rate Q. The cases Q > and Q < 
correspond respectively to injecting or sucking fluid. As 
in the traditional Hele-Shaw problem, the flow in the vis- 
cous fluid is potential, v = V(/i, but now with a velocity 
potential given by [|5| 
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Incompressibility then yields Laplace equation V^(/) — 
for the fleld (p (but not for the pressure). The two bound- 
ary conditions at the interface which complete the defini- 
tion of the moving boundary problem are the usual ones, 
namely, the pressure on the viscous side of the interface 
p = —an, where a and k are respectively surface tension 
and curvature, and the continuity condition for the nor- 
mal velocity u„ = n ■ V(/). The crucial difference from 
the usual case is in the boundary condition satisfied by 
the laplacian fleld on the interface due to the last term 
inEq.(§). 

This problem is well suited to conformal mapping tech- 
niques 0. The basic idea is to find an evolution equa- 
tion for an analytical function z — f{uj,t), which maps 
a reference region in the complex plane uj, in our case 
the unit disk |a;| < 1, into the physical region occupied 
by the fluid in the physical plane z = x + iy, with the 
physical interface being the image of the region bound- 
ary, |a;| — 1. We consider two types of situations, one in 
which the viscous fluid is inside the region enclosed by 
the interface, and one in which it is outside. It can be 
shown [la] that the evolution equation for the mapping 



f{ijj,t) in the rotating case can be written in a compact 
form as 
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where fl* — fe^fi^p/12/i, do = b'^a/12fi, and where 
we have specified the mapping function at the unit 
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The curvature is given by k 



—lm{d'if/{d(j,f\d^f\)} and the Hilbert transform H^ is 
defined by 
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In Eq.(pf), ^ = +1 and § = —1 correspond to the cases 
with the viscous fluid respectively inside or outside the 
interface. 

Explicit time-dependent solutions of Eq. (0) are known 
only for the case d^ — and 57* = 0. Here we report 
explicit solutions for Q.* ^ 0. A class of solutions in this 
context is defined by a functional form of the mapping 
(with a finite number of parameters) which is preserved 
by the time evolution. The problem is then reduced to 
a set of nonlinear ODE's for those parameters. All so- 
lutions we have found fit into the general class of the 
rational form 
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although not any mapping of this form is necessarily a 
solution. A more detailed study will be presented else- 
where jT^]. The general structure Eq.(Q) is known to 
yield explicit solutions also in the nonrotating case (with 
different ODE's for the parameters). However, some im- 
portant classes of solutions of the usual case (fi* — 0) are 
no longer so for fl* ^ 0. This is the case for instance of 
the superposition of a finite number of logarithmic terms 
1^ (sometimes referred to as simple pole solutions, when 
considering the derivative of the mapping rather than the 
mapping itself jl^). These solutions have the interesting 
feature of being free of finite-time singularities 0. On 
the contrary, the 'multiple-pole' case (for the derivative 
of the mapping) turns out to include solutions for both 
il* = and f7* 7^ |0|. Finally, for the polynomial 
case (6j = for all j's), the nonrotating case fi* = 
is known to always yield finite-time singularities in the 
form of cusps [g|,0 . We will see in the rest of this paper 
that this scenario is modified in a nontrivial way by the 
presence of rotation. 

We focus on the role of rotation in preventing cusp 
formation in the subclass of polynomial mappings of the 
form 



j(uj, t) = ao{t)uj^ + an{t)uj 
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For a„ <^ oo this describes a n-fold sinusoidal perturba- 
tion of amplitude a„ superimposed on a circular interface 



of radius ag. It is convenient to introduce the dimension- 
less parameter e = (n-|-z?)a„/ao. The range of physically 
acceptable values of a„ and qq is given by the condition 
< e < 1 for all 71. We also introduce a scaled mode am- 
plitude S — aoan which turns out to be the relevant one 
to characterize the interface instability. To see this, let us 
first compute the standard linear growth rate. Inserting 
Eq.(^ into Eq.(y) and linearizing in a„, we get. 
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The term —Q/2'Kaf^, independent of both n and §, has 
a purely kinematic origin, associated to the global ex- 
pansion (or contraction) of the system. It can easily be 
shown that this quantity would be the growth rate of a 
mode corresponding to a redistribution of area given by 
the undistortcd flow field with radial velocity v = Q/2'Kr. 
This in turn would imply a„(i)ao(t) = const [1^]. Ac- 
cordingly, the marginal modes for 6 (which in the rotating 
case may occur for all n) will be such that the fiow field is 
undistorted by the interface perturbation, although such 
perturbation may grow or decay in the original variables 
a„. Analogously, growth or decay of 5 will correspond 
unambiguously to the stability properties of the actual 
velocity field. In this sense it may be justified to qualify 
the interface instability as described by 5 as 'intrinsic', as 
opposed to the 'morphological' one as described by the 
amplitude a„. In this way the intrinsic growth rate takes 
the simpler form 
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where we have defined Q* = Q /ina^ and have dropped 
the surface tension term, since hereinafter we will restrict 
to the zero surface tension case. We introduce the rel- 
evant dimensionless control parameter of our problem, 
expressing the ratio of centrifugal to viscous forces, as 
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where i? is a characteristic radius of the interface. 

Eq. (0) clearly exhibits the competing effects of rotation 
and injection, although their roles are not quite symmet- 
ric. In fact, notice that Q*, which may have both signs, 
contains a dependence on uq. In practice this means that 
Q* depends effectively on time. An immediate conse- 
quence of this is that the growth of modes is not really 
exponential pq] and may even be nonmonotonic. The 
asymmetry between injection and rotation shows up also 
in the fact that the sign of Q determines which of the 
two effects dominates asymptotically in time. In fact, for 
positive injection rate the typical radius of the inner ffuid 
is growing while typical interface velocities are decreas- 
ing, so centrifugal forces will dominate at long times. On 
the contrary, for negative injection rate, typical veloci- 
ties increase while typical radii decrease, so injection will 
asymptotically dominate over rotation. 



In view of Eq.(^, the most interesting configurations 
will be those in which Q > 0, so that injection and ro- 
tation have counteracting effects. In the case f? = +1 
(viscous fluid inside), which was experimentally studied 
in Ref. Pq] , rotation is always destabilizing. A positive 
injection rate in this case tends to stabilize the circular 
interface. However, for fixed Q, Q* will decrease with 
time, so eventually the interface will reach a radius after 
which all modes are linearly unstable. It is thus expected 
that, in this case, the formation of cusps can only be de- 
layed but not avoided [^ . 

The most interesting case from the point of view of 
preventing cusp formation is i? = —1 and Q > 0, the 
usual configuration in viscous fingering experiments. In 
this case, a small rotation rate will slightly affect the lin- 
ear instability, but will eventually stabilize the growth 
at long times, so it is conceivable to have a nontrivial 
evolution starting from an unstable interface but not de- 
veloping finite-time singularities. 

We now study the fully nonlinear dynamics of polyno- 
mial mappings. Inserting Eq.(||) into Eq.(||) with do = 
we obtain two ordinary differential equations describing 
the evolution of ao(t) and a„(t). These can be integrated 
analytically and yield 
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where kg and A;„ are constants to be determined by initial 
conditions, and where n > 2 for "i? = +1 and ti > 3 for 
i9 = -l. 

Physically acceptable solutions require that the points 
in the cj— plane where d^f{uj,t) — (noninvertible) 
should lie outside the unit disk. The occurrence of a 
cusp is associated to such a point crossing the unit circle 
cjI = 1 at a finite time tc, that is. 
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If we take the initial value ao(0) as the characteristic 
length R, which coincides with the radius of the per- 
turbed circle if we are in the linear regime, and define 
the dimensionless time r = r2*i, condition (O) reads 
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Our aim is now at finding conditions such that an ini- 
tially smooth interface remains smooth for an infinite 
time. Thus we have to impose that Eq.(pT|) should not 
have any solution for Tc > 0. The transition between the 
regions with and without cusps will be defined by the 



conditions that the two members of Eq.(ll) and their 
derivatives with respect to time are equal, such that the 
two curves have a common tangent. These two conditions 
allow us to eliminate Tc, and yield 



X log X — X = — a„fco 
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where x 



2a„ 



and with R 



ao(0) in Eq.(gD. We 
now search for solutions of Eq.(|lJ). For i) = -\-l \t can 
be proven ||l^ that this equation has no solutions, and 
therefore all initial conditions must eventually develop a 
cusp at finite time, as expected from the linear analysis. 
On the other hand, iox ■& — —1 the quantity on the rhs 
of Eq. (|lj) takes the simple form 
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and a nontrivial critical line Pc{e;n) can be found for 
each n > 3. This implies that, in the configuration with 
the viscous fluid outside, for any initial condition there 
is always a certain rotation rate above which there is no 
cusp formation. The numerical determination of these 
curves is shown in Fig.l. 




FIG. 1. Critical lines Pc for different values of n. The re- 
gion free of cusp singularities for a given n is the one above 
the correponding curve. 

The leading behavior for initial conditions in the lin- 
ear regime, e ^ 1, can be found by expanding the Ihs of 
Eq.(|l4|) around x = e and is given by Pc ~ ^^^"^^ ■ No- 
tice that there are qualitative differences for small values 
of n. For n — 3 the curve starts horizontal at the linear 
level, implying that a very small rotation rate is sufficient 
to prevent cusp formation. For n = 4 the threshold curve 
starts with a finite slope and for n > 4 it has an infinite 
slope at e = 0. A more detailed description and analysis 
of this diagram will be presented elsewhere ||l^ . An ex- 
ample of rotation preventing cusp formation is shown in 
Fig.2. 
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sic marginal stability of the circular shape, Q* = 17*. 
Therefore, the sufficient condition, valid for all initial 
conditions of the form Eq.dq), for not developing cusp 
singularities is that a circular interface with radius given 
by ao(0) be intrinsically stable, in the sense of Eq.(R). 
Whether deeper consequences can be drawn in a broader 
context from this inner connection between the linear 
problem and the possibility of cusp formation is an in- 
teresting open question. 
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FIG. 2. Evolution of the interface in the case •& = —1 (vis- 
cous fluid outside), with n = 3, ao(0) = 1.0, e(0) = 0.5, for 
(a) Q.* — 0.025 (cusp formation) and (b) Q,* = 0.045 (cusps 
prevented by rotation). 

In Fig. I we also see that for any given e, the critical 
Pc increases monotonically with n. If we take the limit 
n — > cxD at fixed e we get Ofn^o ~^ 1- From Eq.(14) this 
implies x — 1 and consequently we obtain an absolute 
upper bound p™"^ — 1 for all values of n and e. This 
implies that, for all initial conditions, there is a critical 
rotation rate 
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above which cusps are always eliminated. Although 
Eq.(n6[) has been derived for the class Eq.(p[), it is ex- 
pected that the existence of a certain ilc and the scaling 
with physical parameters given by Eq. dlq) could be more 
general. Notice that P = 1 corresponds to the intrin- 
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